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Passive turbulent drag reduction techniques are of interest as a cost effective means to 
improve air vehicle fuel consumption. In the past, rigid surface waves slanted at an angle 
from the streamwise direction were deemed ineffective to reduce skin friction drag due to 
the pressure drag that they generate. A recent analysis seeking similarities to the spanwise 
shear stress generated by spatial Stokes layers suggested that there may be a range of 
wavelength, amplitude, and orientation in which the wavy surface would reduce turbulent 
drag. The present work explores, by experiments and Direct Numerical Simulations (DNS), 
the effect of swept wavy surfaces on skin friction and pressure drag. Plates with shallow 
and deep wave patterns were rapid-prototyped and tested using a drag balance in the 
Txll inch Low-Speed Wind Tunnel at the NASA LaRC Research Center. The measured 
drag offset between the wavy plates and the reference flat plate is found to be within the 
experimental repeatability limit. Oil vapor flow measurements indicate a mean spanwise 
flow over the deep waves. The turbulent flow in channels with flat walls, swept wavy walls 
and spatial Stokes spanwise velocity forcing was simulated at a friction Reynolds number of 
two hundred. The time-averaged and dynamic turbulent flow characteristics of the three 
channel types are compared. The drag obtained for the channel with shallow waves is 
slightly larger than for the flat channel, within the range of the experiments. In the case 
of the large waves, the simulation over predicts the drag. The shortcomings of the Stokes 
layer analogy model for the estimation of the spanwise shear stress and drag are discussed. 


Nomenclature 


Ax DNS Mesh cell size in the x-direction 

Ay DNS Mesh cell size in the y-direction 

Az DNS Mesh cell size in the z-direction 

Az e Distance along the z-axis between the edge of the sub-domain with undamped sinusoidal wall profile 

and a point located on its outside 
Sq Momentum thickness 

S v Friction length scale i ; fu T 

A Wall wavelength in the direction perpendicular to the wave crests 
Wall wavelength in the direction parallel to the x-axis 
A^ Wall wavelength in the direction parallel to the z-axis 
v Kinematic viscosity 

cr ^3 (i = 1,2,3) z-Components of the stress tensor on the faces of an infinitesimal volume perpendicular to 
the x, y and z axes 

r z z-Component of the wall shear stress 
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var(x, £, T) Localized variance 

£ Coordinate perpendicular to the wall wave crests 

A Curve fit parameter for the drag coefficient (Cd = ARe B ) 

Assl Amplitude of the spanwise velocity oscillations at the wall for the Spatial Stokes layer 
B Curve fit parameter for the drag coefficient (Cd = ARe B ) 

Cd Total drag coefficient 

Cf Friction coefficient 

c p Pressure coefficient 

D Detection function for the local variance 

h Channel half-height 

h w Wall wave amplitude 

k x Wall wavenumber in the direction parallel to the x-axis 

k z Wall wavenumber in the direction parallel to the z-axis 

kth Factor for the threshold of localized variance leading to an event detection 

L x Length of the simulation domain (in the x-direction) 

L z Width of the simulation domain (in the z-direction) 

M Mach number of the freestream 

rii (i = 1, 2, 3) x-, y-, and z-components of the unit vector perpendicular to the wall 
R Ideal specific gas constant for air (287 J/(kgK)) 

Re/m Reynolds number per unit length (^oo/^oo) 

Re^ Reynolds number defined as u^h/v^ 

Re r Friction Reynolds number (u T h/ v) 

Reg Momentum thickness Reynolds number (^ 00 ^ 0 /^ 00 ) 

T Integration time for the calculation of the localized variance 
t Time 

to Detected time of maximum localized variance 

T re f Reference temperature (300.67iL) used in DNS 
u Streamwise velocity component 

u' Streamwise velocity fluctuation with respect to the time-averaged streamwise velocity 

u' rrns Root mean square of the streamwise velocity fluctuations 
Uqq Freestream velocity 

u r ,ref Shear velocity for the reference flat channel 
u T Shear velocity 

v Vertical velocity component 

v' Vertical velocity fluctuation with respect to the time-averaged vertical velocity 
w Spanwise velocity component 

x, y and z Cartesian coordinates oriented in streamwise, vertical and spanwise directions, respectively 
y w Vertical coordinate of the bottom wall of the wavy channel (y w = h w sin(k x x -f k z z)) 


I. Introduction 

T urbulence over air vehicles increases the skin friction drag that in turn drives up energy consumption 
and combustion emissions. In parallel to active flow control and novel aerodynamic designs, passive 
Turbulent Drag Reduction (TDR) methods that rely on incremental infrastructure changes are desirable. 
Riblets and engineered surface roughness address this issue but still face challenges in terms of manufactura- 
bility, durability, and maintenance. Therefore, there is interest in investigating alternate or complementary 
passive techniques to reduce turbulent drag. 

Spanwise traveling waves and surface oscillations were shown to be effective to reduce drag at low Reynolds 
number. 1-5 Viotti 6 showed by Direct Numerical Simulation (DNS) of channel flows at a friction Reynolds 
number Re r = u T h/v = 200 (where u T is the friction velocity of the reference flow, h is the half- height of 
the channel, and v is the kinematic viscocity) that in the near-wall region, the effect of steady streamwise 
oscillations of spanwise velocity forcing in the form of Spatial Stokes Layers (SSL) is similar to that of tem- 
poral oscillations of the wall in the spanwise direction. The net power savings from the SSL were found to 
be superior to the temporal forcing case, owing to an additional spanwise Reynolds stress production term 

2 of 23 


American Institute of Aeronautics and Astronautics 



of negative sign originating from the streamwise gradient of the spanwise velocity. 7 Furthermore, Viotti 6 
showed in selected cases that the scaled DNS phase- averaged turbulent spanwise SSL flow was identical to 
the laminar solution, and used a linear analysis to estimate the net power gain by spatial Stokes layers 
of small characteristic thickness. The practical realization of wall motion poses a number of challenges. 5 
Seeking a passive solution, Chernyshenko 8 proposed to use swept wavy surfaces to generate steady spanwise 
shear stress oscillations. He considered such oscillations to be the dominant drag control element for the 
SSL. He concluded from a linear analysis that a small percentage of total drag reduction is possible for a 
passive swept wavy wall. At Re r = 200, maximum net power reduction was predicted for waves having a 
streamwise wavelength of about 1500 spatial wall units (defined as 5 V = v/u T ), oriented at an angle of ca. 
37° from the mean stream, with a wave amplitude expected to generate an SSL-equivalent peak spanwise 
velocity of the order of 2 velocity wall units (u T ). Swept surface waves were already proposed in the past. 
In Ref. [9,10], Sengupta and Lekoudis proposed a reduced two-dimensional curvilinear turbulent boundary 
layer calculation to model the drag over rigid and traveling swept wall waves. The study indicated a slight 
improvement in skin friction drag with respect to the flat plate in a range of the wave angle, but it was 
over-compensated by the pressure drag. 

The objective of the present work is to revisit the effect of rigid swept waves oriented at an angle from 
the freestream on turbulent drag, and specifically to explore the model proposed in Ref. [8] . To this purpose, 
plates with different patterns of waves were rapid-prototyped and tested using a direct drag balance in a 
low-speed wind tunnel. Three-dimensional Direct Numerical Simulations were conducted on channels with 
flat, wavy and spatial Stokes wall boundaries. Section II describes the wind tunnel tests. Section III provides 
a description of the simulation set-up and selected channel cases. Section IV reports the drag measurements 
and oil vapor flow test results. The simulation results are provided and discussed in Section V. Concluding 
remarks are given in Section VI. 

II. Experimental Set-Up 

A. Tunnel and Drag Balance 

The experimental portion of the study was conducted in the NASA Langley Research Center (LaRC) 7x11 
Inch Low Speed Wind Tunnel (7x11). The 7x11 tunnel is a closed- loop, fan-driven, unpressurized, low- 
speed facility with a top speed of approximately 50 m/s. Broadband turbulence intensity in the freestream 
is approximately 0.5%. Testing was conducted on the lower wall of the rectangular test section that was 
equipped with a single component, floating-element drag balance. The boundary layer on the lower wall of 
the test section was tripped to turbulence with a 0.91 mm diameter rod bonded to the contraction surface 
578 mm upstream of the test model leading edge. Fully turbulent flow was observed for stream velocities 
greater than approximately 7 m/s. The useful unit Reynolds number range for the tunnel is approximately 
0.5 x 10 6 < Re/m < 3.3 x 10 6 . The momentum thickness Reynolds number range at the mid-point of the 
test section is approximately 750 < Ree < 5200. 

The drag balance consists of a linear air bearing, a 0.25 Newton servo-type load cell, and associated 
hardware for leveling and positioning the test surface models. A schematic of the drag balance components, 
as situated beneath the test model, is shown in Fig. 1. The drag balance is fully mechanically isolated from 
the wind tunnel structure. This was primarily required as a means of preventing small, static tilt angles 
caused by tunnel structural distortion that would introduce drag measurement errors. While the drag bal- 
ance is designed to measure horizontal forces, it is also very sensitive to inclination and the streamwise tilt 
angle must be tightly controlled. Tunnel frame distortion can be caused by aerodynamic loading within the 
tunnel loop or by thermally induced deflection of the tunnel structure due to ambient temperature variations 
or tunnel heating. The mechanical isolation also prevents tunnel vibrations from impacting drag balance 
performance. The air bearing system is basically an undamped spring-mass system with a characteristic 
frequency of about 15 Hz. Without the mechanical isolation, the drag balance is subject to large oscillations 
at the resonant frequency rendering the drag data suspect. An earlier version of the system employed a 1/16- 
inch (1.59 mm) thick rubber sheet to bridge the gap between the plenum chamber and support pedestal in 
place of the single pass, non-contact, labyrinth gland as currently shown in Fig. 1. The rubber sheet failed 
to prevent transmission tunnel vibrations and had to be removed. The labyrinth gland currently used has 
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passage clearance of approximately 0.004 inch (0.1 mm). 

The flotation gap around the model was designed as 0.01 inch (0.25 mm). The gap was checked after each 
model change and the model re-centered as required. The test model upper surface was maintained flush with 
the surrounding tunnel wall by use of leveling screws as shown in Fig. 1. The leading- and trailing-edge verti- 
cal misalignments were kept to 0.001 inch (0.025 mm) or better as measured with a dial indicator depth gage. 

Air currents within the flotation gaps and within the plenum chamber are limiting factors in the ultimate 
accuracy and repeatability of the drag balance. The air bearing necessarily discharges air into the plenum 
that must be removed. Also, as the tunnel speed varies, the mean test section pressure drops causing a 
pressure difference between the plenum and test section. This difference causes flow in the flotation gaps 
leading to drag measurement errors. To prevent or minimize gap flow, the static pressure difference between 
the test section and plenum chamber was monitored with a 10 Torr full scale electronic pressure gage and 
adjusted to zero before each drag reading. Pressure equalization was accomplished by an external blower 
at constant rotational speed and a remotely controlled gate-type throttling valve. It was found necessary 
to introduce a small, fixed amount of external ambient air leakage into the plenum to allow the particular 
blower /gate- valve combination used in this study to operate effectively. For very low drag readings where 
the effects of dynamic pressure due to air currents in the plenum are commensurate with test model viscous 
forces, additional design improvements may be required, or, a different technique employed, such as elastic 
support flexures instead of the air bearing. 

The air bearing itself has internal characteristics that can affect its performance. The particular air 
bearing used in the present investigation employs internal discharge orifices located in an offset fashion 
within machined depressions to distribute the air equally between the two exits. This allows the air bearing 
to operate “friction-free” overall but this characteristic is due to large equal and opposite internal viscous 
forces within the small clearance channels of the bearing. Any alteration of the internal equilibrium flow 
distribution not accounted for in the initial tare reading will cause drag measurement errors. Another type 
of air bearing commonly found uses internal porous metals or plastics for air distribution. In that case, 
the airflow pattern may be less predictable and even subject to different start-up flow patterns. The key to 
accurate use of air bearings for viscous drag measurements is to show that errors associated with air bearing 
design and operation are acceptable for the range of model forces being considered. The current design, as 
will be shown in Section IV, has 1-2% repeatability. Improvement beyond this level requires additional study. 

The streamwise test section pressure gradient due to boundary layer displacement thickness growth on 
the four upstream contraction walls and test section walls can also cause flow in the flotation gaps. The 
pressure gradient was approximately adjusted to zero by diverging the rear wall of the test section. An 
additional factor that could conceivably affect repeatability is the assumption that each time the tunnel is 
started, the turbulent boundary layer repeats exactly the same. Given the variability of boundary layer 
transition (only one wall was tripped in the current study), this is a tenuous assumption and may contribute 
to the 1-2% repeatability. Further study of this topic was beyond the scope of the current work. 


B. Test Models 

The flat plate and two wavy wall models were CNC-machined from aluminum using a commercial rapid- 
prototype service provider. Photographs of the models are shown in Fig. 2. The wavy plates were sized 
following the guidance of Ref. [8]. The waves have a peak-to- valley amplitude of 2 h w = 232 /im in one plate 
and 2 h w = 1045 fim in the other. Below, these models will be referred to as the shallow waves and deep 
waves, with notations “WW232” and “WW1045”, respectively. The sweep angle of the waves measured from 
the axial (x) coordinate is 37.4°. The wavelength A in the direction perpendicular to the wave crests is 26.8 
mm for each case, giving aspect ratios of h w /X = 0.43 and 1.95%. The waves are damped within one wave- 
length at the plate periphery, as shown in Fig. 3(b). The net volume of the wavy plates is identical to that 
of the flat plate. The wetted area is increased by 0.013% for the shallow waves and 0.27% for the deeper waves. 

The surface of each model was lightly glass-beaded at the service provider to remove CNC-machining 
marks and subsequently polished at LaRC prior to drag testing. Using a dilute detergent solution as a 
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Figure 1: Schematic representation of floating drag balance components (for illustration only - not to scale). 



(a) Flat plate. 


(b) Wavy plate with h w /X = 0.43%. 


(c) Wavy plate with h w /X = 1.95%. 


Figure 2: Test plates. 
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(a) Plate geometry. 



(b) Normalized surface cutlines through the plate center. 


Figure 3: Plate with swept surface waviness. 


lubricant, the models were first lightly sanded with medium abrasive, conformal foam pads followed by hand 
polishing 600-grit wet carbide sandpaper. The peak-to-valley wave height of each model was sampled with 
a dial-indicating depth gage and found to be within 2% of the design values. 

Data acquisition was accomplished using a digital oscilloscope with 12-bit effective resolution (Agilent 
Model MSO-7034A in high resolution mode). The force sensor output was low-pass filtered at 0.1 Hz, - 
12dB/octave, and unity gain. The sensor output was sampled at 5 Hz and averaged over one 50-second 
oscilloscope trace. The range of forces encountered varied from about 0.02 to 0.2 N on models with a plan 
area of 0.062 m 2 for an average shear stress range of 0.32 to 3.2 Pa. 

Oil vapor flow visualization tests (i.e., the “smoke-wire” method) were performed to check for evidence 
of spanwise deviation of the mean flow over the waves. Prior to the oil test, the plates were painted with flat 
black lacquer to increase the contrast of the oil vapor over the aluminum surface. The oil was lightweight 
mineral oil. The stainless steel (type 304V) smoke wire had a diameter of 0.13 mm and was situated 7 mm 
downstream of the model leading edge, approximately 0.5 mm above the surface. Power was provided by an 
adjustable DC power supply and manually switched through a solid state relay. The vapor was side-lit with 
a cold fiber optic light source. Imaging was accomplished by full frame digital video sequences. 

III. CFD Modeling 

Three-dimensional channel Direct Numerical Simulations (DNS) were conducted using a compressible fifth 
order Weighted Essentially Non-Oscillatory (WENO) code for spatial discretization, with third order Total- 
Variation-Diminishing (TVD) scheme for temporal integration. 11 13 The mesh size is 1224 x 257 x 416 in 
the x, y, and z directions, respectively, with the x-axis parallel to the inflow direction and the y-axis oriented 
in the vertical direction. The mesh is uniform in x and z directions and is refined at the walls using a Roberts’ 
transformation. 14 In the flat channel the cell size along the walls is such that Ax + < 2.8, A y + < 0.18, and 
Az + < 5.8, where the + superscript indicates wall units. At the channel center, A y + < 3.5. 

Five channel cases are studied at Re r = 200, with the following wall boundaries: (1) flat walls with 
non-slip boundary conditions (reference), (2-3) flat walls with spanwise velocity forcing W = AssLsin(k x x ), 
where AssL/^r.ref = 2 and 6, and (4-5) walls with slanted waves defined by y w = h w sin(k x x + k z z)(+2h ) 
with two amplitude-to- wavelength ratios, h w /\ = 0.05 and h w /X = 2.25%. The latter amplitudes corre- 
spond to the SSL velocity forcing amplitudes AssL/u r , r ef = 2 and 9 per the model proposed in Ref. [8]. 
Table 1 summarizes key simulation parameters. The Reynolds number based on the channel half- height 
and freestream velocity is Re^ — 3200. u T , re f is the reference flat plate skin friction velocity that was used 
a-priori to size the wavy wall and Stokes problems. The notation A^ SL refers to the SSL wall velocity 
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Table 1: DNS Parametrization 


Common parameters 

Re r 

200 

Re h 

3200 

M 

0.3 

V* r,ref 

0.062 

Wavy wall 


1520 

kx/^z 

0.7657 

AsSL/u r ,ref 

2,9 

hw / X 

0.5, 2.25% 

X x / h 

7.6 

Lx/Xx-) L z / X z 

2 

Spatial Stokes Layer 


1520 

AsSL/u r , re f 

2,6 

A ssl (actual u T ) 

2.17, 7.65 


oscillation amplitude non-dimensionalized by the SSL skin friction velocity obtained from DNS. The domain 
size is common to all cases and was chosen to comprise two periods of the wavy wall in each direction, as 
shown in Fig. 4(a). There is no phase shift between the top and bottom wall profiles of the wavy channels, 
to avoid mean flow deviations from a varying cross-section. 

Periodic boundary conditions are applied in the stream- and span- wise directions. The flow is maintained 
by applying a body force in the streamwise direction, determined at every step of the time marching by 
requiring that the average mass-flux remains constant. Depending on the wave amplitude, there may be 
a net spanwise component of the mean flow over the wavy wall. In absence of z-directed body force, such 
cross-flow may cause a spanwise pressure gradient affecting the streamwise drag. To assess the effect of 
the boundary conditions, two additional cases were simulated, with only one wavelength of wavy surface 
in spanwise direction, smoothly transitioning to flat edges towards z = 0 and z = L z . Such a channel is 
shown in Fig. 4(b). The transition is achieved by multiplying the channel profile by a cos 2 (k z Az e ) factor 
at its periphery, where Az e is the distance along the z-axis between the edge of the region with undamped 
sinusoidal wall profile and a point located on its outside. 

IV. Experimental results 

Fig. 5 shows the average of two drag data sets for the models WW232, WW1045, and the reference smooth 
flat plate that were acquired over two consecutive days to account for any day-to-day variations in drag. The 
solid curves are least square power laws of the form Cd = AR f where Cd is the total drag coefficient and 
Re is the unit Reynolds number per meter. Curve fit coefficients are listed in Table 1. Fig. 5 shows little, if 
any, significant difference between the reference flat plate and the two wavy wall models. 

The power law curve fit permits closer examination of the drag data relative to the smooth flat plate. 
This is shown in Fig. 6, where the percent difference between the wavy walls and flat plate is plotted based 
on the curve fits in Table 1. The zero value on the vertical axis represents the two-day flat plate average in 
Table 2. The figure shows that the day-to-day non-repeatability for each model’s drag is at about the same 
level as drag differences between the models. For this reason, the small drag differences are not considered 
significant and within this 1-2% band of uncertainty, no definitive statements regarding the drag increase 
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Y 



(a) Channel with wavy surface extending over the 
full xz-domain of area 2X x x 2X Z . 



(b) Channel with center wavy surface of area 2X x X 
1X Z . The waves are damped towards the flat bound- 
aries at z = 0 and z = L z . 


Figure 4: Physical simulation domains for the wavy channels. The black arrow indicates the direction of the 
inflow. 



Unit Reynolds No. (x 10 6 / m) 

Figure 5: Average total drag coefficient with respect to the unit Reynolds number. 
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or decrease of the wavy wall models can be made. The increased uncertainty at lower Reynolds numbers 
is consistent with the discussion of drag balance error sources in Section II. The data show, however, that 
the wavy walls, especially the large amplitude WW1045 model, do not increase the drag beyond this 1-2% 
uncertainty level. 


Table 2: Power law curve fit parameters, Cd = AR^ 


Model 



Day 


A 


B 

Flat plate 



1 


0.109 


-0.247 

55 



2 


0.074 


-0.220 

55 



2- day average 


0.092 


-0.234 

WW232 



1 


0.100 


-0.240 

55 



2 


0.086 


-0.231 

55 



2- day average 


0.093 


-0.236 

WW1045 



1 


0.088 


-0.232 

55 



2 


0.108 


-0.247 

55 



2- day average 


0.098 


-0.240 
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Figure 6: Drag coefficients relative to the smooth flat plat based on power law curve fits. 


A. Oil vapor flow test 

Oil vapor flow visualization tests were performed to check for any spanwise deviation of the mean flow over 
the waves. Fig. 7 shows photographs of oil vapor flow tests. In Fig. 7(a), the amplitude of the flow oscillations 
over the plate with h/X = 0.43% is very small. The mean flow remains mostly two-dimensional, with no 
noticeable spanwise component. In contrast, the streaklines over the plate with h/X = 1.95% in Fig. 7(b) 
show large spanwise oscillations indicative of the non-linearity of the disturbance induced by the waves. 
There is also a net spanwise mean flow, as for example indicated by the deviation of the wave highlighted 
in Fig. 7(c) from the direction of the inflow. The deviation angle is found to vary within 2. 5^1.5°. These 
results support the assumption made at the z-periodic boundaries for the shallow waves. For the deep waves, 
the validity of the periodic conditions used remains questionable. 
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(a) Wavy plate with h w /X = 0.43%. (b) Wavy plate with h w /X = 1.95%. (c) Wavy plate with h w /X = 1.95% (em- 
The red arrow indicates the direction The red arrow indicates the direction of bossed and sharpened version of the image 


of the inflow. the inflow. 


in (b)). The dash-dot line highlights a par- 
ticular streakline. The horizontal dashed 
line shows the position of the wire. The 
three parallel dashed lines indicate the di- 
rection of the inflow. 


Figure 7: Oil vapor flow test. 


V. CFD Analysis 

A. Average turbulent flow 

Fig. 8 shows streamlines of time- averaged velocity at the surface of the bottom wall and in the center hori- 
zontal plane ( y = h/ 2) of the wavy channel with h w /X = 0.05%. A mean flow deviation of less than 0.1° from 
the x-axis (towards negative z) is observed. Although not shown here, the same observation is made for the 
domain with flat edges per Fig. 4(b). This minor deviation is of higher order than the oscillatory flow per- 
turbation, and is not expected to affect the drag analysis. In the case of the large waves with h w /X = 2.25%, 
Fig. 9(a) reveals flow along the wave troughs. At the bottom wall, the convergence of streamlines along 
the downhill slope of the waves, and their divergence back to a locally parallel flow over the uphill slope is 
indicative of flow lift-up and lift-down, respectively. This is illustrated in the streamwise velocity profiles 
of Fig. 11. A flow along the wave troughs is also observed in the wavy area of the domain with flat edges, 
Fig. 10(a). Additionally, Figures 9-10(b) show that there is a net spanwise mean flow component. The 
transition length of X z /2 from waves to flat edges used in the domain of Fig. 4(b) is too short to allow the 
flow to realign along the domain edges. For both partially and fully wavy domains, the angle formed by the 
streamlines with the x-axis at mid-channel height is about 2.6°, which is within the range of the oil vapor 
flow measurements for h w /X = 1.95%. 

The average skin friction coefficient of the flat channel was computed as Cf — 7.46 x 10 -3 , which is in 
reasonable agreement with the empirical expectation Cf = ( 2 i?eI)V 4 = 7.89 x 10 -3 . 15,16 Table 3 summarizes 
the drag results obtained for the different test cases after a minimum of 5 flow passages across the channel, 
normalized by the flat channel data. In this table, a positive drag reduction means that the drag is reduced. 
The drag in the SSL channel with A^ SL = 7.65 is 33.9% lower than for the flat channel, which compares 
well with the TDR of 32% obtained at A,*, = 1250 in Ref. [6]. The SSL channel with A^ SL = 2.17 yields 
8.5% drag reduction. The wavy wall leads to 0.3% skin friction increase and 1.1% pressure drag, yielding 
1.4% total drag increase. This result in within the experimental resolution window of 2% reported in Section 
IV. The skin friction drag increase for the deep waves is found to be 2.7%, with an excessive pressure drag 
that is still being analyzed in the light of the mean flow drift observations. The drag data obtained from the 
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(a) Lower wall surface. 


(b) Channel center (y = h). 


Figure 8: Streamlines of average velocity in the wavy channel with h w /X = 0.05%. The flow goes from left to 
right. 



(a) Lower wall surface. (b) Channel center (y = h). 


Figure 9: Streamlines of average velocity in the wavy channel with h w /\ = 2.25%. The flow goes from left to 
right. 
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(a) Lower wall surface. 


(b) Channel center (y = h). 


Figure 10: Streamlines of average velocity in the wavy channel with h w /X — 2.25% and flat edges at z = 0 
and z — L z . The flow goes from left to right. 



x/h^u/uoc 


Figure 11: Streamwise velocity profiles (solid lines) at different positions along the wavy channel (of which 
the boundaries are shown by dashed lines) with h w /X = 2.25% at z — 0. 
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wavy domains with flat edges at z — 0 and z = L z compare to those for the fully wavy domains. Neither the 
wind tunnel tests nor the DNS models were able to demonstrate turbulent drag reduction for the selected 
wall and flow parameters. 


Table 3: Drag with respect to the flat channel 


Drag component 

SSL A + = 7.65 

SSL A+ = 2.17 

h w /X = 0.05% 

% Friction drag 

66.1 

91.4 

100.3 

% Pressure drag 

- 

- 

1.1 

% Drag reduction 

33.9 

8.6 

-1.4 




Figure 12: Comparison of the surface shear stress computed by DNS at z — \ z to that predicted by linear 
theory, in the spatial Stokes layer (LHS) and in the wavy wall (RHS). The wavy wall profile y w scaled by a 
factor of 1/50 h is shown in the right panel. 


In earlier work, a good agreement was reported between the spanwise wall shear stress simulated by DNS 
in spatial Stokes layers, and that predicted by the linear model presented in Ref. [6,7]. The left panel of 
Fig. 12 compares the spanwise shear stress per unit of spanwise forcing predicted by this model, to that 
obtained by the DNS of the two SSL layers considered here. The agreement on the amplitude is within 
0.5% for A^ SL = 2.17 and 8.5% for Ag SL = 7.65. The linear model proposed in Ref. [8] for the wall shear 
stress was based on a thin boundary layer assumption with a pressure perturbation of the same wave number 
as the wall, and no-slip boundary conditions at the wall. The right panel of Fig. 12 shows a comparison 
between DNS and this model. The vertical gradient of spanwise velocity is predicted within 20% accuracy. 
However, the maximum amplitude of the total z-component of the wall shear stress r+ = Y^i=i (where 

are the z-components of the stress tensor at the wall and rii are the components of the unit vector 

perpendicular to the wall surface) is 50% smaller than on the downhill slope and 25% smaller on the 
uphill slope. This is a result of the non-planarity of the wall and contributions from the streamwise and 
vertical velocity components to the spanwise shear stress. Following the calculations reported in Ref. [8], a 
decrease in spanwise shear stress amplitude by 50% or more with respect to the optimum value would lead to 
less than 1% in net power gain by the wavy wall. In Ref. [9,10], Sengupta and Lekoudis proposed a turbulent 
curvilinear boundary layer model for rigid or moving swept wavy walls, that accounts for the coupling of 
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the wave induced shear stress and pressure gradient. It uses the Loyd model for turbulence closure of the 
disturbance equation system, 17-19 solved consistently with the 2D mean flow equations. The wave geometry 
was similar to that studied here (with 2.54 cm wavelength and h w /X = 0.2%), but the Reynolds number was 
much larger (Reg = 2.25 x 10 6 ). No net drag reduction was found for the rigid waves at any sweep angle. 
The power analysis of Ref. [8] could be revisited with the help of such model validated by DNS simulations, 
to establish whether there is a range of swept wall wave geometries and Reynolds numbers leading to some 
percentage of net drag and power reduction. 





Figure 13: Time- and space-averaged profiles of streamwise velocity u , streamwise velocity fluctuations u rrns , 
Reynolds stress —u'v', and spanwise velocity uJ J. 


In the following, an overline denotes a time-averaged variable (/ = /q^°° f(t)dt) and a prime denotes the 
fluctuation of a quantity around its time-averaged value (/' = f(t) — /). The term “phase-averaged” refers 
to the integration in both time and space over a periodic domain f(y ) = /q^°° dt f Q x dxf Q z f(x,y, z,t)dz. 
u, f, and w are the streamwise, vertical and spanwise velocity components, respectively. Fig. 13 shows 
a comparison between profiles averaged over the full xz-domain (2X x , 2X Z ) in the bottom half (y < h ) of 
the channel with flat wall, wavy walls, and Stokes layers. The streamwise wall shear stress, RMS of the 
streamwise velocity fluctuations u' rrns = and the Reynolds stress component — u'v' are all decreased 

for the SSL. The z-vorticity component for A$ SL = 7.65 is decreased, as expected from the significant 
decrease of the shear stress at the wall. The peaks of u' rrns and shear stress are shifted by about 5 wall units 
away from the wall, as reported in prior work. 6 The phase- averaged velocity profile for the channels with 
wavy walls closely overlap that of the flat channel near the wall, as expected from the previous drag results. 
The maximum of u' rrns at y + — y+ = 15 is 8.5% larger for the shallow waves than for the flat channel. The 
maximum of u' rrns is 8.5% lower for the deep waves than for the flat channel, with a flat maximum extending 
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into the log layer. The Reynolds stress for the deep waves peaks at y+ = 50. Figures 14-15 show turbulence 
statistics relative to the flat channel, along a line perpendicular to the wave crests at y + —y+ = 15 and along 
the wall. The Reynolds stress —u'v' is maximum on the uphill wave slope. No clear trend is seen for u' rrns 
over the shallow waves, whereas maxima and minima are found close to = 0 on the positive and negative 
slopes of the deep waves, respectively. The wall shear stress and pressure coefficients reach their maximum 
value on the uphill slope, approximately three-quarters of a wavelength out of phase with the wall geometry. 




Figure 14: Deviation of u rrns , —u'v' , and Cf in the wavy wall with h w /\ = 0.05%, relative to the flat wall case, 
along a cutline perpendicular to the wave crests. The wave profile y w /h w and pressure coefficient normalized 
by its maximum value c p /c p ,max are also shown in the right panel. 


B. Instantaneous flow 

Fig. 16 shows the streamwise component of vorticity in a yz-plane. The x-vorticity is visibly decreased for 
the Stokes layer and increased for the large waves. Figures 17-18 show 3D top views of a bottom section of 
the channel (y < h, z < A^/4) at constant u' + = (u — u)/u T = —2.0, with time steps of 5 computational 
units, where a computational unit is defined as h/ yjRT re f, with R = 287 J/(kgK) and T re f = 300.67 K. 
The phase-averaged skin friction velocity is used for scaling. Streaks convect in the near wall region at a 
velocity between 8 and 12 wall units (u T ), consistent with Ref. [1]. The SSL shows reduced streak density 
compared to the flat channel. The streaks for the shallow wall waves show similarities to the flat channel in 
terms of their large scale dimension and density, with a higher density of small-scale structures than in the 
flat channel. Such structures can also be seen for the deep waves. One might speculate whether they partly 
result from the convection of streaks through regions of different average velocity profiles while passing over 
wall troughs and crests. 

In order to gain insights into the flow dynamics, we examine next the evolution of instantaneous stream- 
wise velocity profiles around events of peak velocity fluctuation in the buffer zone. Such events are detected 
using the following criterion for the localized variance Far(x, t,T), proposed by Blackwelder and Kaplan in 
Ref. [21] 


where k t h is a threshold parameter, for example 1.2. T is the averaging time for this Variable-Interval- 


1 f t+ ^ / 1 

var(x,t,T) = - j , u 2 (x,s)ds~ I- J , u(x,s)ds 


> k t hu: 


1 2 

rms 


(i) 
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Figure 15: Deviation of u rrns , —u'v', and c/ in the wavy wall with h w /X = 2.25%, relative to the flat wall case, 
along a cutline perpendicular to the wave crests. The wave profile y w /h w and pressure coefficient normalized 
by its maximum value c p /c p ,max are also shown in the right panel. 



(c) h/X = 0.5% (d) h/X = 2.25% 


Figure 16: Instantaneous streamwise vorticity in the plane defined by x — L x /24. 
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(c) to - 10, t+ - 40.4. (d) to - 10, t+ - 26.3. 



(i) to + 5, + 20.2. (j) to + 5, t^ + 3.2. 



(k) to + 10, t+ + 40.4. 



(1) to + 10, t+ + 26.3. 



(m) to + 15, t 60.6. (n) to 15, t 39.5. 


Figure 17: Top view of isocontours of stream wise velocity perturbation u + m —2.0 for the flat channel (LHS) 
and the SSL with Ag SL = 7.65 (RHS). The flow goes from left to right. 
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(a) t 0 - 15, tj 


- 63.8. 


(b) t 0 - 15, t+ -86.18. 



(c) t 0 - 10, t+ - 42.5. 



(e) to — 5, to” — 21.3. 



(k) t 0 + 10, t+ + 42.5. 



(d) t 0 - 10, t+ - 57.45. 



(f) t 0 - 5, t+ - 28.73. 



(j) *o + 5, to" + 3.2. 



(1) t 0 + 10, t+ + 47.45. 



(m) to + 15, t^J” 63.8. (n) to H - 15, t^~ 86.18. 


Figure 18: Top view of isocontours of stream wise velocity perturbation u'+ — —2.0 for the wavy wall with 
h w / X = 0.05% (LHS) and the wavy wall with h w /X = 2.25% (RHS). The flow goes from left to right. 
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Time- Average (VITA) calculation. We use here the same value T + = 10 as in Ref. [21], where again the 

subscript indicates wall units (=v/v%. for time). In the limit of T -A oo, Far(x,£, T) = u® ms . When an 
event is detected per (1), the detection function D{t + ) = 1.0. It is equal to zero otherwise. The detection 
coordinate is y + = 15 for all cases, except for the large Stokes waves for which we use y + = 20, according 
to the 5 wall unit shift of the peak of u' rrns observed in Fig. 13. It is worth noting that the phase- average 
of u^ ms was used for detection, instead of the (x, z)-dependent time average. This allows avoiding multiple 
detections from a low local average value, as can arise for the wavy wall (see Figures 14-15). 

Fig. 19 shows the temporal evolution of the streamwise velocity perturbation u' and localized variance 
uar(x, t + ,T + ) at sampled locations within the domain of Figures 17-18, where events with D = 1 were 
detected. The reference time for each case was set to the time of peak variance mr(x,t,T). It is also 
the reference time used in Figures 17-18. A spacing between events of the order of a few hundred wall units 
is consistent with the experimental results of Ref. [21]. The streamwise velocity fluctuation amplitude for 
the spatial Stokes layer with spanwise velocity forcing Ag SL = 7.65 is noticeably lower than in the other 
cases. No clear difference is identified between the flat channel and the shallow waves. The deep waves show 
multiple peaks of localized variance close to the detection threshold. Figures 20-23 show the instantaneous 
and average streamwise velocity for each channel, at times close to the peak of the localized variance. Ideally, 
such profiles would result from averaging over multiple events, using conditional averaging as described in 
Ref. [21]. Given the computational space allocation and time requirements for such an analysis, only in- 
stantaneous profiles are plotted here, bearing in mind that they represent isolated snap-shots. The detected 
events show an inflectional transition from a streamwise momentum defect to a momentum excess, followed 
by an increase of the wall shear stress, and subsequently a recovery toward the average velocity profile. Once 
again, the SSL shows reduced fluctuations around the average profile and no obvious difference is observed 
between the flat channel and the shallow wave cases. The large waves show a large increase in wall shear 
stress at + 5 and somewhat irregular profile slope changes. 


VI. Conclusions 

In this work, the turbulent drag over wall waves slanted at an angle from the streamwise direction were 
studied by experiments and Direct Numerical Simulations. Test plates were rapid-prototyped with two 
wave amplitude-to-wavelength ratios of 0.05% and 2.25%, at a wave crest angle of 37.4° with respect to the 
streamwise direction. They were measured in the NASA LaRC 7x11 Inch Low Speed Wind Tunnel (7x11) 
at a unit Reynolds number 0.6 x 10 6 < Re/m < 2 x 10 6 . The drag force measurement resolution for the 
set-up was estimated as +/-2%. Within this measurement accuracy, no drag improvement or degradation 
could be determined on these models relative to the flat plate. Oil vapor flow measurements showed that the 
flow over the waves with h/X = 0.05% is parallel to the inflow, whereas the mean flow over the deep waves 
shows large oscillations with a mean spanwise component at an angle of about 2.5° from the inflow. 

The turbulent flow in channels with flat walls, swept wavy walls and spatial Stokes spanwise velocity 
forcing was simulated at a friction Reynolds number of two hundred by DNS. The phase- averaged and 
instantaneous turbulent flow characteristics of the shallow wavy wall with 0.05% amplitude-to-wavelength 
aspect ratio were found to compare to those of the flat channel near the wall. The friction drag is maximum 
on the uphill slope of the waves and minimum on the downhill slope. The pressure coefficient is out of 
phase with the wave profile by a quarter of a wavelength. The total drag in the channel with shallow waves 
was computed as 1.4% larger than in the flat channel, within the repeatability limit of the measurements. 
Simulation results were presented for the deep waves, however the assumed spanwise periodic boundary 
conditions need further consideration based on the observed mean spanwise flow. The DNS results for the 
Spatial Stokes layer were consistent with prior work in terms of drag reduction, turbulence statistics and 
linear model predictions. The linear model based on the spatial Stokes layer analogy predicted the amplitude 
of the vertical gradient of the spanwise velocity over the wavy wall within 20%, but it was found that the 
total spanwise shear stress is at least twice smaller, so that the model does not provide an accurate estimate 
of the net power gain. Existing linear models accounting for the wave induced shear stress and mean flow 
modification could be used in combination with DNS to enhance prior net drag and power calculations over 
an extended range of swept wave geometries. 
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t+ 

(a) Flat channel at the location x = 0.75A;e, z = 0.12X Z . 



t+ 


(b) Wavy wall with h w / A = 0.05% at the location x = 
0.92A*, z = 0.12A Z . 






(c) Wavy wall with h w /X = 2.25% at the location x = (d) SSL with A^ SL = 7.65 at the location x = 0.75A X , z = 


0.5Acc, z = 0.12A Z . 


0.12A*. 


Figure 19: Time evolution of the velocity and localized variance fluctuations, along with the detection threshold 
and function. 
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var(T ), k th u'£ a 



Figure 20: Streamwise velocity in the flat channel at the location (x,z) = (0.75A X , 0.12A^). Red line and 
symbols: instantaneous value; black line: time average. 



Figure 21: Streamwise velocity in the channel with wavy wall with h w /\ — 0.05% at the location (x,z) — 
(0.92ACC, 0.12A^). Red line and symbols: instantaneous value; black line: time average. 
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Figure 22: Streamwise velocity in the channel with wavy wall with h w /X = 2.25% at the location (x,z) = 
(0.5Acc, 0.12A Z ). Red line and symbols: instantaneous value; black line: time average. 



Figure 23: Streamwise velocity in the SSL channel with Ag SL = 7.65, at the location (x, z ) = (0.75A X , 0.12A Z ). 
Red line and symbols: instantaneous value; black line: time average. 
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